Transition voltage spectroscopy: a challenge for vacuum tunneling models at nanoscale 



O 
(N 

oo 






c^ 



loan BaldcEO and Horst Koppcl 
Theoretische Chemie, Universitdt Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany 

Several recent studies on the transition voltage (Vt) of molecular and vacuum nano-junctions are 
based on calculations of the tunneling current through an energy barrier within the so-called Sim- 
mons model, which is an approximate WKB-type approach developed for thin insulating films with 
infinite transverse extension. In this paper devoted to vacuum nano-junctions, we compare the Sim- 
mons results for Vt with those obtained from the exact Schrodinger equation by exactly including the 
classical (non-retarded) charge image effects. The comparison reveals that the Simmons estimates 
for Vt are completely unacceptable for nanogap sizes (d) at which image effects are important. The 
Simmons treatment drastically overestimates these effects, because it misses the famous 1/2 factor 
related to the fact that the image interaction energy is a self energy. The maximum of the Simmons 
curve Vt vs. d turns out to be merely an artefact of an inappropriate approximation. Unlike the 
Simmons approach, the "exact" WKB method yields results, which qualitatively agree with the 
exact ones; quantitative differences are important, demonstrating that the transmission prefactor 
has a significant impact on Vt- Further, we show that a difference between the work functions of 
(say,) the left and right electrodes, which gives rise to a Volta intrinsic field, may be important for 
the ubiquitous asymmetry of the measured 7-V^-characteristics [I{V) 7^ —I{ — y)] in general, and for 
the different Vt-values at positive and negative biases reported in vacuum nano-junctions in partic- 
ular. The weak dependence Vt = Vt(d) found experimentally in vacuum nano-junctions contrasts to 
the pronounced dependence obtained by including the exact electrostatic image contribution into 
the exact Schrodinger equation. This demonstrates that not only the molecular transport, but also 
the transport through a vacuum nano-gap represents a nontrivial problem, which requires further 
refinements, e. g, a realistic description of contacts' geometry, retardation effects due to the finite 
tunneling time, local phonons, surface plasmons or electron image states. 

PACS numbers: 73.63. Rt, 85.35.Gv 85.65. -l-h. 
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I. INTRODUCTION 

In the continuous efforts for miniaturization, using sin- 
gle molecules as active components for future nanoelec- 
tronic devices appears at present as the only conceivable 
alternative, which escapes the fundamental limitations 
of complementary metal-oxide semiconductor (CMOS) 
technologies. In molecular devices, electron transfer be- 
tween the source {S) and drain (D) electrodes across 
the nanogap (width d) can occur via through-bond and 
through-space processes. In these processes, electrons 
have to tunnel through an energy barrier given by the en- 
ergy offset SB = min {ep ~ £homo,£lumo ~ £f) of the 
closest molecular orbital (HOMO or LUMO) from elec- 
trodes' Fermi level (ep) and the metallic work function 
W, respectively. Similar to electron transfer in numerous 
chemical reactions, through-bond and through-space pro- 
cesses can compete, because (i) in usual off-resonance sit- 
uations, the Fermi level lies close to middle of the HOMO- 
LUMO gap (charge neutrality), and Eb amounts several 
electronvolts, i. e., only slightly smaller than W and (ii) 
charge image effects, which narrow, round, and lower 
the barrierfii^ arc more pronounced in vacuum than in 
molecules due to the different dielectric constants [k,. = 1 
versus k^ ~ 2 — 3, respectively, cf. Eq. ^ below]. Con- 
cerns that the measured currents are not (or not only) 
due to the active molecules have plagued the field of 
molecular electronics since its inception.— For a proper in- 
terpretation of transport experiments, through-bond and 
through-space mechanisms should be each correctly un- 



derstood. 

Transition voltage spectroscopy (TVS) has been pro- 
posed recently as an appealing tool to deduce Sb, a 
key parameter for molecular junctions4r— TVS relies 
upon simple intuitive considerations inspired by the 
barrier picture. The initial TVS conjecture was that 
the minimum y = Vt of the Fowler-Nordheim (FN) 
curves [log{I/V'^) versus l/V, obtained by recasting the 
source-drain /-^-characteristics] suffices to determine 
the molecular energy offset: eVt = £b- This minimum 
was ascribed to the point where the shape of the energy 
barrier tilted by the applied voltage V {eb -^ SB^eVx/d, 
where x specifies the position) changes from trapezoidal 
to triangular. Calculations of the tunneling current 
within the Simmons' model^ challenged the validity of 
the barrier description for molecular junctions,— claiming 
that it yields a dependence of Vt on eb and d disagreeing 
with cxperiments4i"— lis Based on them, it was suggested 
that the dependence Vt = Vt{d^ or the magnitude of Vr^ 
can be used to discriminate between through-bond and 
vacuum tunneling. Some aspects of the TVS in molecular 
junctions^ii^rii and gated single-molecule transistors^^ 
received theoretical consideration. 

Letting alone the fact that the calculations of Ref. [3 
used the Simmons' approximation^ (an issue to be ad- 
dressed in detail in the present paper), it is worth em- 
phasizing that those calculations disregarded an impor- 
tant fact. Namely, that in order to phenomenologically 
model the tunneling through a molecule, it is too sim- 
plistic that, irrespective of the molecule under consider- 



ation, electrons within the barrier ("molecule") be gen- 
erally characterized by the free electron mass m. Several 
studies drew attention on the fact that this is an illegit- 
imate assumptiouji^i^ and one rather needs to employ 
an effective mass m* , which is molecule specific [with the 
present notations m* = 7Tt*(rf, £b)]. For the above rea- 
sons, the utilization of the barrier picture for molecular 
transport is only possible if the phcnomcnological param- 
eters £b and m* are determined by fitting experimental 
data or from companion ab initio calculations. 

While modeling the active molecule as an energy bar- 
rier may be too simplistic, and describing the electrons 
that tunnel in molecules (unless they possess delocal- 
ized TT electrons) within the effective mass approxima- 
tion questionable, the barrier picture appears as the most 
natural framework for studying the transport in vacuum 
nano junctions. Still, it is only recently that a theoretical 
study becomes attractive, because, in contrast to the im- 
mense number of studies on molecul ar j unctions (a very 
incomplete list includes Refs. l3-[7l ll(]|ll8l - l23l and citations 
therein) experiments on electron tunneling through vac- 
uum in nanojunctions have been carried out only very 
recentlyjii 

Among very numerous approaches to the elec- 
tron transport by tunneling relying upon the bar- 



rier picture? 
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an approach due to SimmonsS be- 



came increasingly popular for the interpretation of 
important experimental findings in molecular-/iiaiio- 
j unctions 4i^—iiS'i^iii'i2ii^ This popularity among exper- 
imentalists should certainly be related to the fact that, 
unlike in other (more accurate, see, e. g., Ref. Isil) treat- 
ments, Simmons gave approximate but simple analytic 
formulas, which can be easily used for fitting experimen- 
tal data. While primarily discussing and formulating the 
barrier approach to nanojunctions in general, we will also 
to draw attention on the fact that Simmons' formulas, 
aiming to describe thin insulating films, are not appro- 
priate for molccular/nano-transport. An important issue 
to be emphasized is that certain predictions based on 
these formulas are artefacts, which preclude an adequate 
interpretation of valuable experimental findings. 

The remaining part of this paper is organized as fol- 
lows. In Section |lTl the Simmons approach will be re- 
viewed, emphasizing its limitation. The confinement of 
electron motion in transverse directions, an aspect which 
is important for nanojunctions with atomic contacts but 
is not accounted for within the Simmons approach, will 
be considered in Section IIIII Section |IV] is devoted to 
the comparison of the results for Vt obtained by solv- 
ing the Schrodinger equation exactly with those of the 
"exact" WKB and ( "approximate" WKB-typc) Simmons 
approach, which clearly demonstrates that the last ap- 
proach is inappropriate to account for image effects at 
sizes of interest for nanotransport. In Section |Vl the 
case of electrodes with different work functions will be 
analyzed. An important technical issue, the relationship 
between the apparent barrier height and the electrodes' 
work function, will be addressed in Section I VII Based 



on the present exact results for Vt , in Section I VIII it will 
shown that the validity of the WKB approximation has 
been too optimistically assessed in previous studies. The 
main results of the present study will be summarized in 
Section FVlIII where possible sources for the disagreement 
between the present theory and experiment will be indi- 
cated. 



II. REVIEWING THE SIMMONS' APPROACH 

The theoretical framework of the single-particle tun- 
neling through an energy barrier was established long ago 
by Sommcrfeld.— The current density due to the (elastic) 
tunneling along the x-direction across a sandwich con- 
sisting of an insulating film placed between two (source 
and drain) electrodes caused by a bias V^ > between 
source (chosen below as the negative electrode) and drain 
(positive electrode) with chemical potentials ^5 = and 
ji£) = —eV can be expressed as2 



J 



'/l3 



eV 

eVl T{E^,V)dE, 



£r(p.,F)d3p = — — (1) 
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E,T{E^,V)dE, 
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where T is the transmission coefficient. To get Eq. ([2|), 
we assumed a parabolic conduction band dispersion Ep = 



^{pI+pI+pD-Sf- 



To deduce the current, a very popular approach is to 
employ the WKB expression of the transmission 



TwRsiE^^V) = e 



Jd, 



dx 's/(pBix) — Ej: 



(3) 



Like many others (e. g, Refs. HHHillillOlll) , 
Simmons^ did not use the above "exact" WKB approxi- 
mation, but rather an approximate version of the "exact" 
WKB expression 



Tr{E,,V) = exp 



-A^/eB - Ea 



(4) 



replacing the general position dependent barrier 0b (x) 
by a rectangular barrier. In Eq. ([4]), the fact that 
the barrier 4>Bix) is not rectangular is accounted for 
by considering an effective (average) rectangular barrier 
of height eb and a spatial extension As that can be 
smaller than the distance between electrodes (embedded 
molecule/nanogap size) d as well as a correction factor 7. 
This factor enters the quantity A = 2^As{2my^^/h. For 
certain insulating thin films, Simmons estimated 7 ~ 1 
within at most a few perccnts. A correction factor 7=1 
has been implicitly assumed in Refs. IMllI although sev- 
eral earlier studies drew attention on the fact that this 
assumption is inconsistent with experimental I-V-data 
for molecular transport. Alternatively, in view of the 
above expression of A, departures from a rectangular- 
shaped barrie r^i^i^^ can be considered by employing an 
effective mass m* 7^ m (m — >■ m* = 771,7^) J^i^^ 



By using Eq. (|4]), Eq. ([2]) can be integrated exactly in 
closed analytical form. The resulting, rather lengthy ex- 
pression was given previously by various authors (e. g., 
Refs. I3l26[) including Simmons;^ who only retained the 
leading terms in his numerical calculations. The central 
result of the Simmons approach^ is the following expres- 
sion of the current density J 



J/Jn = ^b{i + 



l-fO 



[a {eb + eV) 



,-Ai 
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O e 



-AieB+ep) 



- {es + eV) 
(5) 



where Jq = e/h/{2n^/^sy'. Eq. ^ applies if the Fermi 
level of (say,) the right electrode lies above the bottom 
of the left electrode, < eV < ep', otherwise the second 
term in the RHS ("backward" current) should be omit- 
ted. This y-range suffices for practical purposes, since it 
is much broader than the range of the biases, which the 
molecular junctions fabricated so far can withstand or 
those used in experiments on vacuum nanoj unctions pii 
Simmons' leading order expression ([5]) was utilized in 
a series of more recent works^"— li^ii^iiiiiSi^ to interpret 
a series of valuable experimental findings. An impor- 
tant pragmatical advantage of the approaches based on 
Eqs. (m and ([5]) is that it can also simply account for 
image effects, sec Eq. (|8])^i^i2^ 

By supposing identical electrodes (work functions 
Ws = Wd = W) first, the total barrier (J)b{x) through 
which electrons have to tunnel can be expressed as^i^i^i 



(I)b{x) = Sb — eVx/d + (j)i{x). 



(6) 



In addition to a bare (rectangular) energy barrier sb 
[eb — W for vacuum tunneling) , the energy barrier com- 
prises the contributions of the applied bias V and of the 
charge images 0i(a;)i^ 

The effective barrier width As = S2 — si is determined 
by Simmons from the barrier extension (0 < Si < x < 
S2 < d) at the electrodes' Fermi energy without bias (no- 
tice that energies are measured throughout relative to 
the Fermi level) 



0s(a;)U= 



0. 



The effective barrier height is expressed as 

1"^ dx 



SB ^ £b - eV 



Sl +52 

2d 



Ai 



(t)i{x). 



(7) 



(8) 



In the range of experimental interest {e\V\ < Eb) and 
ignoring image effects {(f>i = 0), the barrier is trapezoidal, 
Sl = 0, As = S2 = d, and Eb — £b ~ eV/2. 

For infinite planar electrodes placed at x = and x ~ 
d, the electrostatic interaction energy (j)i{x) between an 
electron located at x and its images can be expressed 
exactlji^ 

c2 



(|)^{x) 



4Krd 



-2^(1) -I- V 



+ ^ 1 



(9) 



where ip is the digamma function. To better under- 
stand the later analysis, we include throughout a di- 
electric constant Kr although for vacuum (nr = 1) it is 
superfluous. In order to work out Eqs. ([7]) and ([5|) in 
closed analytic forms, various parabolic approximations 
of the RHS of Eq. ([9|) around its maximum at x = d/2, 
(j), (f ) := -^ log2 were employed in Refs. ilHii. In- 
stead of Eq. ©, Simmons used the "approximate" ex- 
pression 



^hx) 



-1.151og2 
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d- 



which he deduced from his "exact" expression 



^n^) 
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[ndY — xP- nd 



(10) 



(11) 



In fact, Simmons' "exact" expression (|lip is incorrect. 
The series entering the RHS of Eq. (fTT|) can be summed 
out, and the result is twice the RHS of Eq. ©• An easy 
way to understand this wrong factor is to consider sit- 
uations very close to one electrode, e. g., x > 0; then, 
Eq. (fTTI) yields (\)i(x) ~ —e^ji^Krx). This is twice the 
result well known from electrostatics textbooks; the at- 
tractive force between a point charge e placed at x > 
and an infinite plane at origin is F^ix) = — e^/(4KrX^), 
and the corresponding interaction energy is 0i(x) = 
— / d^Fa;(^) = — e^/(4Krx). Simmons' exact expres- 
sion misses nothing but the "famous" 1/2 factor related 
to the fact that the image interaction energy (/)i(x) is a 
self energy and not the charge-potential product)^ 

Recently, Simmons' expression (jlOp has been taken 
over in studies on image effects on the transition 
voltage.i^^ In those works, all the terms (i. e., not only 
the leading ones written explicitly) in Eq. ([5]) were used 
in numerical calculations. To obtain the results pre- 
sented below referred to as the Simmons results, we have 
also used these full expressions, as previously done in 
Refs. [gflll. The terms ©(...) which are not explicitly 
written in Eq. ([S]) may yield "corrections" to Vt up to 
^ 10%f2i, for smaller values of d and 'eb (^s^ ^ 1)- 
Still, we emphasize that the inclusion of these terms is 
a priori questionable, as Eq. ([5]) is basically the result 

of a (simplified) WKB approximation, and a valid WKB 

1 1I 
treatment requires sufficiently large values of Ae-^ (see 

Ref. [Sy and the discussion below). 

As a curiosity, we note that in order to improve the 
agreement between the theory based on the Simmons ap- 
proach and the experimental data for vacuum tunneling, 
a correction factor C, ~ 0.4 — 0.7 multiplying the RHS of 
Eq. (jTOl) has been empirically introduced in Ref. [ill. It 
would be tempting to identify C with the aforementioned 
factor 1/2 missing in Eqs. pT|) and ([TU]). Unfortunately, 
the remedy of the Simmons results is not so simple, as it 
will be shown below. 

Let us summarize the basic features of the Simmons' 
approach: 
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FIG. 1: Energy barrier 0s (x) for eb = W = 5.7 eV, d = 8 A, 
and Kr = 1 computed using the exact contribution of charge 
images, Eq. ^, and using the Simmons' "approximation", 
Eq. HOI), at F = Vt = 2.12 V. Notice that the height and 
width of the Simmons barrier are substantially smaller than 
the exact ones, and therefore image effects are drastically 
overestimated by the Simmons model. 



(i) Electrons tunnel across an energy barrier 0b (cc) 
in x-direction and move freely in transverse (y, z)- 
directions, cf. Eq. ([T]). 

(ii) The transmission coefficient across the barrier is 
computed within an approximation more restrictive than 
the WKB-approximation, which assumes that the barrier 
profile 4>b{x) can be replaced by a rectangular effective 
barrier 



As 



\/4>b{^ 



Ex 



7 



Ves — Ex ~ y/sB — E^, (12) 



where the (Ex- and 1^-independent) correction factor can 
be taken 7 ~ 1. This amounts to consider that, irrespec- 
tive of their energy (fijj < E^ < fJ-s), all electrons that 
contribute to the current tunnel across a barrier of the 
same width As = S2 — si and height eb, which is the 
barrier "seen" by the electrons at the Fermi level of the 
unbiased junction {E.^ = 0), cf. Eqs. ([7]) and ([8]), and 
Fig.ffl 

(iii) The contribution of charge images to the energy 
barrier is estimated from Eq. (1101) . This overestimates 
the image effects by a factor ~ 2.^ 

To demonstrate the fact that the Simmons approach 
is inadequate for nanojunctions, we will compare below 
the results of this approach both with the results of the 
"exact" WKB approximation [i. c., by using Eq. (|3]) and 
not Eq. (|3])] as well as with the exact results obtained 
by solving the Schrodingcr equation with the exact bar- 
rier potential, Eq. ([6]). The exact transmission T(£'x, V) 
can be easily obtained by imposing the overall continu- 
ity of the wave function and its derivative. Inside the 
junction (0 < a; < d), the wave function ^{x) obeys the 
Schrodinger equation 



ft2 ^2 



*(a;) = 0. 



(13) 



The exact Schrodinger equation written above can easily 
solved analytically in the absence of images ((/f^ = 0)i^ 
The exact wave function ^{x) for this case can be ex- 
pressed as a linear combination of the Airy functions 
Ai(z) and Bi{z), where z = ^ — x/\ and 
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The transmission is given by T{Ex , V) 
J\f{Ex,V)/V{Ex,V), where^ 
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d? 



{A'sB'j, - A'^B's) + ksko {AsBd - AnBg) 



^ {AsB'^ A'^Bs) + ^ {AdB's ^sBd) 



Above, the obvious dependence on Ex and V has been 
omitted for brevity, the subscripts S and D correspond 
X = and x = d, respectively, and the prime stands for 
the derivative of the Airy functions with respect to z. 



e. g.. Ad = Ai{C-d/X) and B'g = [dBi{z)/dz]^^^, ks = 

yj2m{£F + Ex)/h,, and kr, = ^2m,{£p + Ex+ eV)/h. 

In the presence of image effects, we have performed 
exact numerical calculations for the potential given by 
Eqs. ([6]) and (|9]) or (fTO)) . To this aim, we used a suf- 
ficiently large number N of small pieces (until reach- 
ing the convergence) = xo < xi < . . . < xn = d 
and considered piecewise constant potentials 4'B,n = 



III. LATERAL CONFINEMENT 

Eq. ([5]) represented the starting point for interpreting 
many valuable experimental results for molecular junc- 
tions. Therefore, it is noteworthy that, in fact, the de- 
scription based on Eq. ([iJ, from which Eq. ([5]) is de- 
duced, is appropriate neither for molecular junctions nor 
for vacuum nanojunctions with atomic contacts, because 
in those cases the electron motion transverse to the junc- 
tion is confined. Rather than being continuous variables 
(0 < Eyz < £p at zero temperature), the transverse 
energies are quantized, and only the lowest-energy trans- 
verse channel contributes to electric conduction. There- 
fore, to describe electric transport in such junctions, the 
integration over the transverse (y, z) directions should 
be omitted, and the counterpart of Eqs. ([T]) and is a 
current expressed by 



'-^-l 



f r{Ex,V)dEx. 

J-eV 



(14) 



It is worth emphasizing that Eq. (|14p qualitatively differs 
from Eq. ^. To understand that Eq. ([2]) is basically in- 
appropriate for nanotransport, let us consider the linear 



conductance 



G ex lini 



dJ{V) 
dV 



TiE,,V^O)dE, 



(15) 



Eq. (J15p is inconsistent with Landauer's fundamental 
statement for the transport at nanoscale ("conductance 
is transmission", G ex T{eF,V ~ 0)). The electrons 
contributing to Eq. (|15p can have a kinetic energy E^ 
of motion across the junction spanning the whole con- 
duction band (0 < E^ < ep); E^ is not restricted to 
Ex ~ ep- This is the consequence of the fact that the lat- 
eral (y,z) confinement is missing in Eq. ([2]). The states 
contributing to linear response in Eq. ([T]) have a total 
energy Ep = E^ + Ey + E^ — ep-, but the transverse 
kinetic energies are continuous variables in the range 
< Ey^z < sp- To conclude, Eq. ^ may be appro- 
priate for traditional electronics envisaged by Simmons, 
but not for nanojunctions with atomic contacts, wherein 
the lateral confinement is essential and Eq. ([Ti)) should 
be used. 

Mathematically speaking, the only difference between 
the situation without lateral confinement discussed in 
Section IIVI and that with lateral confinement is the 
manner of carrying out the energy integration, namely, 
Eq. ([2]) versus Eq. (fT4| . Whether computed exactly by 
solving the Schrodinger equation ([T5|) or approximately 
via Eq. ([3]) or Eq. ^ the transmission coefficient is given 
by the same formula. This is the main reason why the 
shortcomings of the Simmons results presented below 
are essentially the same, irrespective whether the one- 
dimensional or three-dimensional description, underlying 
Eqs. dH]) or ^, respectively is employed. 



Let us first examine the case of a relatively wide 
(d = 8 A) vacuum junction (k^ = 1) with electrodes char- 
acterized by identical work functions Eb = Ws = Wd = 
5.7 eV (i. e., the highest experimental value, which is the 
most favorable for a valid WKB treatment). The consid- 
eration of vacuum nanojunctions obviates an important 
conceptual issue for molecular junctions. In the latter, in 
view of the finite tunneling time (t) spent by electrons 
within the barrier ("molecule"), it is unclear whether the 
static dielectric constant Hr{^)\u]=o or the dynamical di- 
electric Kr{'jj)\u)'^i/T should be used in Eq. ^. 

Results for the vacuum nanogap specified above are 
presented in Fig. [2l The lower group of three curves 
(which can hardly be distinguished among themselves 
within the drawing accuracy) represent exact, WKB 
and Simmons results obtained by ignoring image ef- 
fects. For the bias range of experimental interest {V < 



/(nA) 




2.5 v{y) 



IV. SIMMONS' RESULTS VERSUS WKB AND 
EXACT RESULTS 

In the cases typical for thin dielectric/semiconducting 
films, with (practically) infinite transverse extension and 
widths d ~ 20 — 50 A, large dielectric constants k^ ~ 10 
[cf. Eqs. ^ and (|10l) ]. and voltages that are not too high 
analyzed by SimmonSji^ the approximations underlying 
Eq. (O may be reasonable. 

However, we will show that the above approximations 
are not justified for molecular/nano-junctions. Although 
we are going to present only results obtained by account- 
ing for the lateral confinement discussed in Section lllll we 
note that there are only insignificant quantitative differ- 
ence between the exact. WKB. and Simmons results also 
when the lateral confinement is ignored. To emphasize 
again, the shortcomings of the Simmons approach dis- 
cussed in this paper are not specific for the cases where 
the lateral confinement is important. 

To obtain the numerical results presented below, we 
fixed the Fermi energy to the value ep = 5.5 eV spe- 
cific for gold, employed experimental work function val- 
ues Ws,D = 3.2 — 5.7 eV (see also Section |Vl] below), and 
considered nanogap sizes d relevant for cxperimentSf^ 



FIG. 2: /-^-characteristics without and with images. Notice 
that tlie curve deduced from full quantum-mechanical cal- 
culations (label "exact") without and with images has been 
multiplied with a factor 0.250 and 0.503, respectively, which 
represents the deviation of the ohmic conductance of the "ex- 
act" WKB value from the "exact" conductance. 

3V)^ ignoring image effects, the Simmons approxima- 
tion agrees well with the "exact" WKB approximation. 
Concerning the latter, it basically deviates from the ex- 
act quantum-mechanical calculations by a constant fac- 
tor lexact — -IbQIw K B ■ Although this constant fac- 
tor precludes a quantitative analysis of experimental I- 
y-characteristics, it does not notably affect the transi- 
tion voltage {Vt = 2.13 V for the lower group of curves 
without images in Fig. [5]), deduced from the minimum 
of log(//l/2)(= log/ - 2 logy). For the highest work 
function Ws — Wd — 5.7 cV deduced in expcrimentsjii 
a difference 5Vt ~ O.IV (comparable to the experimen- 
tal inaccurac y"*^' ^ ^ i ^-^ ) from the approximate estimates and 
the exact value sets the smallest nanogap size dmm — 8 A 
for reliable WKB and Simmons approximations for this 
case. Briefly, the Simmons (and WKB) approximation 
yields reasonable estimates of the transition voltage Vt for 
barriers sufficiently high and wide if image effects were 



negligible. However, this situation deteriorates at smaller 
sizes (see Fig. [5^) and heights even without charge im- 
ages. 
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FIG. 3: Exact, WKB, and Simmons results on the transition 
voltage for vacuum tunneling (a) without and (b) with charge 
images. In the latter case, results obtained within the Sim- 
mons approach are shown both with the incorrect and with 
the correct image forces [Eqs. (|10|l and (|9]), respectively]. The 
curves for r = Aeb^ /4: reveal that the condition r > 1 of 
Ref. [3^ for a valid WKB approximation is too optimistic. 

The situation drastically changes when image effects 
are considered. Because they rounds off the barrier cor- 
ners, diminish and narrows the barrier hill;^ image effects 
increase the tunneling current. Image effects are most 
pronounced in vacuum (lowest k,, ~ 1). For the case of 
Fig. [21 the exact current is enhanced by three orders of 
magnitude. Due to the incorrect factor in Eq. (|lip im- 
age effects are drastically exaggerated by the Simmons 
approach. The example of Fig. [2] demonstrates that not 
only the Simmons results, but also the WKB values sig- 
nificantly deviate from the exact ones. They become 
unacceptable at the relatively large barrier widths and 
heights at which they would be reasonable without image 
effects. This holds true not only for the magnitudes of the 
currents, but also for the values estimated for the transi- 
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FIG. 4: Transition voltage for vacuum tunneling without im- 
ages computed exactly, and within the WKB and Simmons 
methods for electrodes with equal and different work func- 
tions Ws,D specified in the legend. Without image effects, 
the Simmons curves for {Ws = 5.7 eV, Wd ~ 3.2 eV) and 
{Ws = 3.2 eV, Wd = 5.7eV) coincide. Notice that, with 
our choice, the source/drain (S/D) is the negative/positive 
electrode and Vt-values are always positive. 



tion voltage. This is illustrated by the exact, WKB and 
Simmons results for the transition voltage in the pres- 



encc of charge images, which are presented in Figs. [SJd 
and m As visible there, Simmons' estimates for Vt are 
acceptable only for large sizes {d > 20 A) where image 
effects are negligible altogether. 

Still, this is not the whole issue concerning the va- 
lidity of the Simmons model in cases where image ef- 
fects are significant. Most importantly, the Simmons ap- 
proach even yields a qualitatively incorrect prediction. 
This fact becomes clear by inspecting the curves depicted 
in Fig. [5}d, and by comparing the panels of Fig. [5] among 
themselves. As visible there, the Simmons curves are 
qualitatively incorrect; they exhibit a maximum, which 
the exact curves do not display. The trend toward 
smaller sizes is just opposite: while, roughly, Vt is in- 
versely proportional to d at large sizes in both cases, at 
low sizes the exact Vt-values increase superlinearly with 
l/d whereas the Simmons estimates decrease with 1/d. 
Simmons curves exhibi ting a maximum have been previ- 
ously shown in Refs. lolllll. where the three-dimensional 
description [Eq. ^] was implicitly assumed. The present 
comparison with the exact curves demonstrates that this 
maximum is nothing but an artefact of the Simmons 
approximation, and this applies irrespective whether a 
lateral confinement is accounted for or not. The occur- 
rence of this Simmons maximum for curves of Vt{d) is 
an artefact of the WKB rectangular-barrier approxima- 
tion [approximation (ii) of Section |TT], while the fact that 
this maximum is located at such small sizes is the con- 
sequence of the incorrect factor of Eqs. ((TT|) and ^TU\\ (as 
noted under (iii) in Section |Tl|) . Computations based on 
the approximate transmission given by Eq. (j4]) using the 
image potential energy of Eq. ^ instead of Eq. ([TU)) still 
yield a maximum in Vt (d) albeit located at substantially 
smaller sizes. This fact is illustrated by Figs. I^b, and 
[5t,d. 



V. ASYMMETRY DRIVEN BY DIFFERENT 
ELECTRODES' WORK FUNCTIONS 

So far, we have considered electrodes characterized by 
identical work functions {Ws — Wd)- This is a reason- 
able assumption for vacuum tunneling between planar 
crystalline electrodes with macroscopic transverse exten- 
sions consisting of the same metal (even) if their separa- 
tion (d) falls in the nanometer range. However, this is 
an unlikely situation in molecular/nano-electronics with 
atomic contacts. Even for electrodes of the same chemi- 
cal nature (e. g., gold), it is unlikely that the local crystal 
orientation at contacts, which is hard to control experi- 
mentally, is identical. This general feature has been con- 
firmed by the experimental data from the vacuum tun- 
neling junctions fabricated in Ref. I 111 , for which values 
VF ~ 3.2 — 5.7 eV have been deduced (see also Section 
IVII below) . This fact strongly supports the idea that not 
only the work functions characterizing different junctions 
are different, but also the values characterizing the (say,) 
left and right electrodes of a given junction can be dif- 



ferent. This effect seems to be not properly considered 
in recent studies on nanoelectronic devices with atomic 
contacts, but may be quite relevant, as already pointed 
out by Sommerfeld [see Ch. 20(b) and Fig. 36 of Ref.Q. 
In the case of electrodes with different work func- 
tions {Ws 7^ Wd), a nonvanishing Volta potential dif- 
ference arises, which yields an intrinsic (Volta) field 
{Wd — Ws)/{ed) that favors the motion of electrons from 
the metal with the higher work function to that with 
the lower work function)^ In such situations, the above 
Eq. © should be generalized to express the total barrier 
(J)b{x) to account for the Volta field as foUows^i^i^i 

0s (.t) =eB + {Wd ~Ws- eV) x/d + 0,(x). (16) 

It is worth emphasizing that the Volta field is significant: 
the above example reveals that values of Wd — Ws can 
be a few electronvolts, such that the Volta field can be 
even larger than the applied field, since molecular/nano- 
junctions can hardly withstand voltages larger than V ^ 

More or less asymmetric experimental I-V- 
characteristics for positive and negative biases [i. e., 
^{^) 7^ ~^(~^)] ^I'S ubiquitous in the transport at 
nanoscale. This is also the case of the vacuum nano- 
j unctions of Ref. [ll|. A consequence of this symmetry 
breaking is the fact that the transition voltages for 
positive and negative biases V are different, a fact 
clearly visible in Fig. 1 of Ref. [ll|. 

One can easily convince oneself that the I-V- 
characteristics deduced from Eqs. (]§]), (HH) [or ([T|)] by 
using the barrier © are symmetric, I{V) = —I{—V). As 
a consequence, the transition voltages deduced for pos- 
itive and negative biases are of equal magnitude. How- 
ever, this symmetry becomes broken in the cases where 
the work functions are different. As an illustration of 
this asymmetry, in Figs. |4] and [5l we present results ob- 
tained by employing the experimental values of the work 
functions {Ws,^^^)} = {5.7, 3.2} eV mentioned above, 
both by considering and ignoring image effects. The 
curves for the transition voltages computed exactly ex- 
hibit the aforementioned asymmetry both without and 
with charge images, in qualitative agreement with exper- 
iment. 

The work function asymmetry {Ws 7^ Wd) can repre- 
sent a general source for asymmetric /-V^-characteristics 
for nano junctions and even for devices based on sym- 
metric molecules. In molecular devices, it can further 
enhance the asymmetry of the electric potential profile 
obtained by solving the Poisson-Schrodinger equations 
selfconsistently. 

As in the case of identical electrodes, the Simmons 
curves for Vt are unacceptable. By ignoring image effects, 
the Simmons curves are not asymmetric (cf. Fig. HJ;), 
while by considering image effects all these curves exhibit 
a maximum (cf. Fig. [SJ;) , which is absent for the exact 
curves and is nothing but an artefact of an inadequate 
approximation. 



VI. APPARENT BARRIER HEIGHT VERSUS 
WORK FUNCTION 

In the calculations presented above we have employed 
the (extreme) values of the work functions Ws,d (3.2 eV 
and 5.7 eV) given in Ref. ^^ These values have been de- 
duced there from tunneling conductance G measurements 
by assuming a dependence G ex exp(— 2(iv2rnA/?i) at 
larger d^^ and identifying the apparent barrier height A 
with W. In Fig.ini we present results for the conductance 
computed exactly with and without charge images. In 
the absence of image effects, the curves G{d) demonstrate 
a virtually perfect exponential decay down to d — >■ and 
the deduced apparent height is A = W = {Ws + Wd) /2. 
By including the image forces, an exponential decay is 
still visible at larger d. However, at smaller d-values 
deviations from the exponential decay are clearly visi- 
ble, which set in the faster, the lower the work func- 
tions are. The apparent height deduced from the large-d 
portion is no more equal to the average work function, 
5 A = A — W y^ 0. From our exact numerical data, we 
found a difference ^A, which does not depend on Ws,d 
and amounts SA ~ — 0.57eV. To further check that this 
difference is a genuine image effect, we also performed 
similar calculations by artificially increasing the dielec- 
tric constant Kr from the vacuum value (k^ = !)• In 
agreement with the presence of Kr in the denominator 
of the RHS of Eq. ([9]), we deduced from our numerical 
results a virtually perfect scaling SA{Kr) ~ — 0.57 eV/Kr- 

Our exact numerical results indicate that even in the 
presence of image forces, the work function asymmetry 
SW = Ws — Wd 7^ has practically no effect on the 
linear tunneling conductance; it solely depends on the 
average value W. In contrast to the case at high voltages, 
the most exotic asymmetric values of Ws.d can hardly 
change the linear conductance by one percent, so we can 
safely ignore this effect at low voltages V. Still, this 
behavior indicates that the linear conductance data do 
not suffice, and a further refinement is needed to estimate 
the values of Ws,d from the experimental data. 

Attempting to undertake this effort seems to make 
little sense at present. Although relatively large, the 
above <5A-value is still smaller than the experimental in- 
accuracy of ^ 1 eV of the apparent barrier heights de- 
duced in Ref. 11. Therefore, in our numerical calcula- 
tions, we have simply taken over the experimental values 
Ws,D = 3.2; 5.7 eV deduced without applying corrections 
due to the fact that 5A^Q and 6W ^0. 



VII. ON THE VALIDITY OF THE WKB 
METHOD 



The results for the transition voltage presented above 
showed that, qualitatively, the "exact" WKB curves for 
Vt behave similar to the exact ones; quantitatively, dif- 
ferences at smaller sizes are significant. Concerning the 
WKB method, wc make the following comment at this 



point. The exact counterpart of Eq. Q, the transmis- 
sion through a rectangular barrier hill (eb > E^) is well 
known 



T{E^,V) 



ikskoKg / K^iks + ko) cosh^(K_Bs) 



{ksko - k|) sinh (kbs] 



(17) 



where kb = [2m{eB — Ex)Y^^/h and ks,D have been 
defined above. [Notice that for rectangular barriers, the 
Simmons and WKB expressions coincide, cf. Eqs. ^ and 
(|4]).] Based on the fact that coshz « sinhz sa exp(z)/2 
for large z, it was claimed that the WKB results should 
not notably differ from the exact ones i^ 



A 



'/4 > 1- 



(18) 



However, by inspecting the curves for r in Fig. [3l one 
must conclude that Eq. (fT8|) does not guarantee reli- 
able WKB and Simmons estimates for Vt- Although 
the hyperbolic functions entering Eq. (J17p reduce to 
an exponential for large arguments, there still remains 
an important difference between the exact and WKB 
transmission, namely the preexponential factor, which 
is energy- and F-dependent. Were this factor constant, 
IwKB = const X lexact, the WKB and exact FN-curves 
log^wKB/V^) = logi^Iexact/V^) -f const would havc the 
minima at the same (transition-)voltage, but our results 
of Figs. |3l m and [5] contradict that y^^^ represents 
a reliable estimates for y^^°-'^* , So, the folkloristic dic- 
tum on a slowly varying prefactor multiplying a rapidly 
varying exponential does not hold for cases of interest in 
molecular/nano-transport. The need to reconsider the 
prefactor in the FN-tunneling theory for field-emission, 
a topic of traditional vacuum electronics, has also been 
pointed reccntly»^ As is well known from textbooks, the 
validity of the WKB approximation is difficult to assess 
in general. Conclusive evidence only comes from the di- 
rect comparison with the exact results, largely rendering 
the WKB calculations superfluous. 

Although the WKB curves for Vt can significantly de- 
viate from the exact ones, at least they are qualitatively 
correct; they do not exhibit the spurious maximum of 
the Simmons curves. Therefore, one may still ask what 
is the key difference between the the WKB and Simmons 
approximations. In the latter, the actual barrier (/)_b(x) 
is replaced by a rectangular barrier, and this comprises 
two aspects. A rectangular barrier does not only mean 
a constant height sb [as if the potential drops at the 
contacts were the same (namely, V/2)], but also a con- 
stant (energy-independent) barrier width s ~ S2 ~ si, 
cf. Eq. 0. Let us examine Fig. [TJ Within the WKB 
method, depending on their energy, the electrons with en- 
ergies in the Fermi window {fis = > Ex > fJ-D = —eV) 
that contribute to the tunneling current "see" a bar- 
rier (f>Bix) — Ex whose width X2{Ex) — xi{Ex) is Ex- 
dependent. Here, xi^2 are defined by 0s(a;i_2) = Ex- 
Within the Simmons approximation, electrons tunnel 
across a rectangular barrier eb — Ex whose width As = 



S2 — Si is the smallest; Xi{Ex) < si and X2{Ex) > S2- 
Practically, all electrons that contribute to the WKB- 
current tunnel through a barrier broader than the Sim- 
mons barrier, and this is the most important reason why 
the Simmons method (even letting alone the missing fac- 
tor 1/2) overestimates the image effects. The approxi- 
mation X2{Ex) — xi{Ex) ~ S2 — Si may be justifiable for 
low biases {fis — fj-o = eV ^ Eb), but at higher biases 
the width of the effective barrier is significantly larger 
than S2 — Si- As a remedy, one can determine the effec- 
tive width by using the barrier midway between fxg and 
Hd, i- e., using (l>Bix)\x=Xi,2 = (ms + M-d)/2 = -eV/2 
instead of Eq. (O. This attempt is a trade off of en- 
hanced (reduced) tunnel probabilities of electron states 
below (above) the orange line of Fig. [U which wipes out 
the spurious maximum of the Simmons curves discussed 
above, but it still represents an approximate treatment 
more restrictive than the "exact" WKB approximation 
ofEq. ©. 



VIII. SUMMARY AND OUTLOOK 

The results reported in the present paper can be sum- 
marized as follows: 

(i) We have drawn attention on the fact that the uti- 
lization of the barrier picture for molecular or nano- 
junctions with atomic contacts should properly account 
for the confinement of the transverse electron motion, a 
feature not included into the Simmons formula ([5]). 

(ii) We have presented a detailed comparison of the 
results for the transition voltage obtained exactly and 
within the Simmons model. Particular emphasis has 
been laid on the fact that the Simmons method, whose 
formulas have been utilized to process valuable experi- 
mental data in several recent studies, is inadequate for 
molecular/nano-junctions. Letting alone that the Sim- 
mons model (a) ignores the transverse confinement, (b) it 
drastically overestimates the image effects because of the 
missing factor 1/2 in the employed expression of the im- 
age interaction energy, and (c) predicts a spurious max- 
imum in the curve for Vt{d), which is an artefact of a 
specific rectangular barrier approximation. 

(iii) We have presented results demonstrating that, as 
far as the transition voltage is concerned, the validity of 
the WKB method was too optimistically assessed previ- 
ously. 

(iv) We have discussed that a nonvanishing difference 
bewteen the electrodes' work functions, which is likely 
the case in most experimental setups based on atomic 
contacts, may be significant for the ubiquitous asym- 
metry of the /-^-characteristics and the different Vt- 
magnitudes at opposite bias polarities. 

Despite intensive and extensive theoretical and exper- 
imental efforts in the last decade, the field of molecular 
transport is still confronted with a series of difficulties. 
The transport problem is alleviated when electrons tun- 
nel through a vacuum nanogap, and one needs no more to 



consider, e. g., how the positions and widths of the molec- 
ular levels change upon contacting to electrodes, what is 
the role played by molecular correlations, what is the po- 
tential profile across the embedded molecule, etc, and all 
these in a nonequilibrium current-carrying state. There- 
fore, aiming to gain further insight into the transport at 
nanoscale by (also) studying vacuum nano junctions, as 
recently done experimentallyfii appears as a noteworthy 
attempt. 

Although, for reasons like those enumerated above, the 
transport through vacuum nanogaps should be easier to 
understand, it turns out to be far from trivial. The 
present study demonstrates that the curves of V* = Vt (d) 
reported for the vacuum nanojunctions investigated ex- 
perimentally in Ref. [ll| are rather challenging. In the 
present paper, we have demonstrated that the initial ten- 
tative attempt to explain the weak and broad maximum 
(almost a broad plateau) of the experimental Vt (d)-curves 
within Simmons-type calculationsii (letting alone that 
the Simmons maxima are much more pronounced and 
narrow than the experimental ones) cannot be justified; 
the computed maxima turn out to be an artefact of the 
Simmons approach. While being able to indicate the 
drawbacks of the Simmons approach, we could not ex- 
plain the weak d-dependence of the experimental curve 
for Vt , although our results for Vt (d) have been obtained 
by solving the Schrodinger equation (|13p exactly and em- 
ploying the expression ([9]) of the charge image contri- 
bution, an exact result of classical electrostatics. The 
weakly d-dependent experimental values of Vt, which fall 
in range ^ 1.5 — 2.5 V (depending on the apparent bar- 
rier height) for d < 8Aii are significantly smaller than 
the theoretical estimates of Fig. [S^. This clearly indi- 
cates that other effects, not included within our treat- 
ment, should be considered, and in the following we will 
enumerate several possibilities. 

Since the "exact" Eq. ^ applies to the case of in- 
finite planar electrodes, described as classical contin- 
uous media rather than atomic arrays, one could at- 
tempt to employ a more realistic geometry, or to con- 
sider effective image planes at a; = xq and x = d — xq 
(with xq ^ 2 A), displaced from the nominal surfaces at 
X = and a; = d, as suggested from intuitive reasons^ 
or by LDA-calculations4i The classical result expressed 
by Eq. ^ can be deduced microscopically as an effect 
of surface plasmons in metallic electrodes4^ By consid- 
ering the static limit (w -^ 0) of the surface plasmon 
polarization, one obtains deviations from Eq. ^ 



(t)i{x) 



4:Krd 






Bu[l--^,0 



2^(1) 



21og(l-w) 



where u = exp(— 2fccd), fee being a surface plasmon mo- 
mentum cutoff, and B is the incomplete beta function.— 
Eq. ^ can be obtained from the above equation by tak- 
ing the ("local") limit /c^ ^- oo. Retardation effects 
(w 7^ 0) due to local phonons, surface plasmons^i^ or 
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finite tunneling timo^ could also play an important role. 
Last but not least, in view of the well known fact that 
the image potential can bind electrons near the surface 
of liquid helium)^ surface image electron states should 
also deserve consideration. A possible role of electron 
states localized at electrode's surface has already been 
suggested in Ref. [llj and discussed more quantitatively 
in Ref. M. 



Acknowledgment 



The financial support for this work provided by the 
Deutsche Forschungsgemeinschaft (DFG) is gratefully ac- 
knowledged. 



* Also at National Institute for Lasers, Plasma, and Radi- 
ation Physics, ISS, POB MG-23, RO 077125, Bucharest, 
Romania 

^ A. Sommerfeld, Z. Phys. 47, 1 (1928). 

^ A. Sommerfeld and H. Bethe, Handbuch der Physik, 
edited by Geiger and Scheel, volume 24, page 446, Julius- 
Springer- Verlag, 1933. 

^ A. Aviram and M. A. Ratner, Chem. Phys. Lett. 29, 277 
(1974). 

* J. M. Beebe, B. Kim, J. W. Gadzuk, C. D. Frisbie, and 
J. G. Kushmerick, Phys. Rev. Lett. 97, 026801 (2006). 

^ S. Ho Choi, B. Kim, and C. D. Frisbie, Science 320, 1482 

(2008). 
^ J. M. Beebe, B. Kim, C. D. Frisbie, and J. G. Kushmerick, 

ACS Nano 2, 827 (2008). 
^ B. Kim, J. M. Beebe, Y. Jun, X.-Y. Zhu, and C. D. Frisbie, 

J. Amer. Chem. Soc. 128, 4970 (2006). 
® J. G. Simmons, J. Appl. Phys. 34, 1793 (1963). 
^ E. H. Huisman, C. M. Guedon, B. J. van Wees, and S. J. 

van der Molen, Nano Letters 9, 3909 (2009). 
^° S. H. Choi et al., J. Amer Chem. Soc. 132, 4358 (2010). 
" M. L. Trouwborst et al., Nano Letters 11, 614 (2011). 
^2 M. Araidai and M. Tsukada, Phys. Rev. B 81, 235114 

(2010). 
" J. Chen, T. Markussen, and K. S. Thygesen, Phys. Rev. 

B 82, 121412 (2010). 
^^ T. Markussen, J. Chen, and K. S. Thygesen, Phys. Rev. 

B 83, 155407 (2011). 
^^ L Baldea, Chem. Phys. 377, 15 (2010). 
^^ R. E. Holmlin et al., J. Amer. Chem. Soc. 123, 5075 

(2001). 
^"^ W. Wang, T. Lee, and M. A. Reed, J. Phys. Chem. B 108, 

18398 (2004). 
^^ P.-W. Chiu and S. Roth, Appl. Phys. Lett. 92, 042107 

(2008). 
^^ H. Song et al.. Nature 462, 1039 (2009). 
^° H. Song, Y. Kim, H. Jeong, M. A. Reed, and T. Lee, J. 

Phys. Chem. C 114, 20431 (2010). 
^^ A. J. Bergren, R. L. McCreery, S. R. Stoyanov, S. Gusarov, 

and A. Kovalenko, J. Phys. Chem. C 114, 15806 (2010). 
^^ H. Song, M. A. Reed, and T. Lee, Adv. Mater. 23, 1583 

(2011). 
^'^ M. C. Lennartz, N. Atodiresei, V. Caciuc, and 

S. Karthaeuser, J. Phys. Chem. C , doi: 

10.1021/jp204240n. 



^* R. H. Fowler and L. Nordheim, Proc. Roy. Soc. London A 
119, 173 (1928). 

^^ L. W. Nordheim, Proc. Roy Soc. London A 121, 626 
(1928). 

^® R. Holm, J. Appl. Phys. 22, 569 (1951). 

^■^ J. Bardeen, Phys. Rev. Lett. 6, 57 (1961). 

^^ W. A. Harrison, Phys. Rev. 123, 85 (1961). 

^^ R. Stratton, J. Phys. Chem. Solids 23, 1177 (1962). 

'■^^ T. E. Hartman, J. Appl. Phys. 35, 3283 (1964). 

^^ K. H. Gundlach, Solid-State Electronics 9, 949 (1966). 

^^ W. F. Brinkman, R. C. Dynes, and J. M. Rowell, J. Appl. 
Phys. 41, 1915 (1970). 

^^ V. B. Engelkes, J. M. Beebe, and C. D. Frisbie, J. Amer. 
Chem. Soc. 126, 14287 (2004). 

^^ J. G. Simmons, J. Appl. Phys. 34, 2581 (1963). 

^^ A. Hartstein, Z. A. Weinberg, and D. J. DiMaria, Phys. 
Rev. B 25, 7174 (1982). 

^® K. Gundlach and J. Simmons, Thin Solid Films 4, 61 
(1969). 

^^ The factor 1.15 entering Eq. (|10p represents an empirical 
correction, which makes the "approximate" average value 
closer to the "exact" one. Without this factor, the max- 
imum values obtained from Eq. (jlip and Eq. (|10|) coin- 
cide [(jii{x = d/2) = — 2e^ Iog2/(«;rd)], but near electrodes 
(i > and X < d) the "approximate" values given by 
the latter equation are smaller than the "exact" values 
given by the former equation. A significantly better ap- 
proximation than Eq. (|10l) . which recovers the exact re- 
sults in the limits a:: > 0, x- ~ d/2, and x < d is^ 



U-) - -^ (4^^^^ + I + I^ 



■^^ Notice, however, that the Volta field only affects the energy 

barrier and not the integration limits in Eqs. ([T]) and (|14|l . 
^^ G. Binnig, N. Garcia, H. Rohrer, J. M. Soler, and F. Flores, 

Phys. Rev. B 30, 4816 (1984). 
■"^ R. G. Forbes, J. Appl. Phys. 103, 114911 (2008). 
''^ N. D. Lang and W. Kohn, Phys. Rev. B 7, 3541 (1973). 
^'^ Z. Lenac and M. Sunjic, Nuovo Cimento B 33, 681 (1976). 
'''^ A. Tagliacozzo and E. Tosatti, Physica Scripta 38, 301 

(1988). 
"'* A. Hartstein and Z. A. Weinberg, J. Phys. C: Solid St. 

Phys. 11, L469 (1978). 
"^ C. C. Grimes, T. R. Brown, M. L. Burns, and C. L. Zipfel, 

Phys. Rev. B 13, 140 (1976). 



11 



V,{V) 
5- 



_ W^W, =i.2eV 



exact w. images 




0.1 



0.15 



0.2 



l/d{A) 



3- 



2-_ 



- lV^=J.2eV;W^=5.7eV 



WKB w. images 




l/d{K) 



V(V) 
3- 



Simmons w. images, Eq. (10) 



- Wj,^J.2eV; W^:^5Je\ 

- Wj,^5.7eV; W^=i.2eV 
_ W =W =i.2eV 




0.1 



0.15 



0.2 



l/d(K) 



V(V)r 

3^ 



1- 



- H'^ = W,=5,7eV 




i/rf(A) 



FIG. 5: Transition voltage for vacuum tunneling consider- 
ing image effects computed exactly, and within the WKB and 
Simmons methods for electrodes with equal and different work 
functions Ws,d specified in the legend. Simmons curves com- 
puted both with the incorrect [Eq. (jlOp . panel c] and with 
the correct [Eq. ((9]), panel d] image forces are shown. Notice 
that, with our choice, the source/drain (S/D) is the nega- 
tive/positive electrode and Vt-values are always positive. 
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FIG. 6: Normalized conductance G/Go (Go = 2e^ /h) com- 
puted with and without image effects versus nanogap width 
d for several apparent barrier heights A given in the legend. 
Without images A = W , while with images A ~ H^ — 0.57eV. 



